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Abstract 

We investigate the perturbation in the tunneling current caused by non- 
adiabatic mechanical motion in a mesoscopic tunnel junction. A theory in- 
troduced by Caroli et al.EHU is used to evaluate second order self-energy cor- 
rections for this non-equilibrium situation lacking translational invariance. 
Inelastic signatures of the mechanical degrees of freedom are found in the 
current-voltage I(V) characteristics. These give rise to sharp features in the 
derivative spectrum, d 2 I/dV 2 . 

PACS numbers: 73.40.Gk 72.10.Di 85.30.Mn 85.42. +m 
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I. INTRODUCTION 



Electron tunneling has provided valuable insight into the physical properties of solid state 
systems. Measurement of the current through stationary tunnel junctions has been used to 
elucidate the density of states in superconductors and many-body properties of metals and 
semiconductors. Through much theoretical effort, a satisfactory theoretical description of 
the tunneling process in condensed matter has been attained. Recently, microfabrication 
techniques have progressed to the point where it is now possible to make extremely compliant 
vacuum tunneling electrodes which may not remain mechanically stationary during the 
tunneling process. In the following work we consider the effect that such motion within 
an elastic tunneling barrier may have upon electron tunneling characteristics. Since, to 
our knowledge, previous studies of this question have been phenomenological@~i, we shall 
attempt to present a more complete theoretical approach to this problem. 

We investigate the characteristics of the tunneling current through a square potential 
barrier where the barrier is fixed on one side and is allowed to oscillate freely on the other. 
The goal of this work is thus to predict the effect of a mechanically compliant electrode, 
which can recoil from a tunneling process, on the current- volt age characteristics and its first 
and second derivatives. Mathematically this results in the treatment of a localized phonon 
representing the movable part of the barrier. 

In order to treat this problem in a many body approach which can be extended to 
more realistic situations than the simplified model considered here, we use a theory treating 
translationally non-invariant systems under non-equilibrium situations developed in a series 
of papers by Caroli, Combescot, Nozieres and Saint- JamesEHl (hereafter referred to as CCNS) 
based on the Keldysh non-equilibrium perturbation formalism0. The theory gives a rigorous 
derivation of an energy dependent transfer term from first principles, thereby extending the 
first, more phenomenological approach of Bardeeni, and it allows for a treatment of the 
phonon perturbation in the usual diagrammatic theory. 

The principle introduced in Ref. [I] can be used to treat a non-equilibrium situation such 



2 



as occurs in metal insulator metal (M-I-M) tunneling. It consists of making one or more 
partitions in the system, allowing for a separate treatment of regions in equilibrium. The 
separate parts are then joined through an appropriate transfer term. For an arbitrary one- 
electron potential Caroli et al. showed that this method yields an exact treatment. CCNS 
then derived in Ref. ^ the well-known expression for the tunneling current /(V) through a 
rigid square barrier for a quasi-equilibrium situation. Even though the series of papers starts 
with a discrete model, CCNS later make a transition to a continuous representation which we 
use for the calculations in the present paper. They then consider electron-phonon interaction 
effects in two following papersi'i, but the discussion of the effects remained incomplete since 
a treatment of real phonons is rather involved. The reason is that, for many body potentials, 
the evaluation has to be done more carefully, as renormalizations due to many body effects 
will always depend on the whole system rather than just pieces of it. However in special 
cases such as the one presented here, it can be shown that a simple approximation will work 
very well. 



II. OUTLINE OF THE CALCULATION 

Our physical system consists of a mechanically-compliant cantilevered metal tip of mass 
m c , placed a small distance 2a from a stationary bulk metal counter-electrode (see Fig. 1). 
The movable tip assembly, which we shall refer to as the "cantilever" , is modeled as a spring 
with Hooke's law force constant k c . Electrical contacts are made to the tip and the counter- 
electrode. A voltage V is applied across these electrodes and the resulting characteristic of 
the current I is measured as a function of the applied voltage. In our model Hamiltonian 
of the system, the metal tip is considered to be a single-mode harmonic oscillator with 
characteristic frequency u c . Assuming a quasi equilibrium situation of equal Fermi energies 
(e r F ~ e l F ) and the simple case of equal work functions for both sides, the resulting potential 
barrier can be approximated by a square barrier with barrier height V$. The oscillation of 
the cantilever simply adds an additional degree of freedom, resulting in a modulation of the 
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barrier width on one side (see Fig. 2). Hence our quantum mechanical Hamiltonian reads 

H =^ + ^+ 2 kcX * + Vo ^ Xe + a) ~ d{ ~ Xe + Xc ~ ° )] (1) 

where the subscripts e and c refer to the tunnelling electron and the cantilever respectively, 
while x c is the instantaneous displacement of the cantilever from its equilibrium position, 
and p c is the corresponding momentum. Also 9 is the usual unit step function. 

To treat the additional term in perturbation theory the potential term is expanded to 
first order in x c , which yields V r (x e ,x c ) = Vo[9(x e + a) — 9{x e — a)] — x c Vo5(x e — a). In 
order to apply a many body treatment to this problem we follow the approach of Ref. |2] and 
split the system into two halves by adding an infinite potential barrier in the middle of the 
unperturbed barrier. We thus obtain two sub-parts similar to the appendix of Ref. ||| with 
the potential terms V r (x e , x c ) = V 9(—x e + a) — x c V 5(x e — a) for x > and V r (x e , x c ) = oo, 
for x < for the right hand side medium and Vi(x e ,x c ) = Vq9(x £ + a) for x < and 
Vi(x e ,x c ) = oo, for x > for the left one (see Fig. 3). As in Ref. || we shall assume that 
these two resulting sub-parts are semi-infinite in the direction perpendicular to the barrier, 
while at the same time having a finite cross sectional area on which we impose periodic 
boundary conditions. In second quantized notation the Hamiltonian for the right side reads 
k 2 - 1 

Hr = Y,77^4 C <l + U; ^ a+a+ o) ~ X oV C i,9|| C 9 ',«||0S.( a )09 '( fl )( a+ + a ) ( 2 ) 

q /me A 1x,q xl ,q\\ 

where the Xq = ^Jl/2m c cu c , the root mean square zero point displacement of the cantilever, 
and we use h = 1. Also the variable x labels the coordinate perpendicular to the barrier. In 
this representation the wave functions <p qx are the x-dependent components of wave functions 
Xq which diagonalize the system when the unperturbed barrier is included; we will discuss 
their explicit form later. We will follow the procedure used in Ref. ^| to obtain an expression 
for the current across the barrier I(V), see their Eq. (41), (46), (47), and for dl/dV , 
d 2 I/dV 2 as well as corrections to these quantities, following the procedure in Ref. ^. As 
indicated in our model Hamiltonian, we assume that the barrier width does not fluctuate 
spatially — this leads to perfectly specular transmission. The unperturbed Green's functions 
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for both parts of the system have been calculated in the appendix of Ref. || and will be used 
as the starting point for further calculations. 

Following Ref. ^, to calculate the current we first obtain a self energy correction for 
the unperturbed Green's functions of our system, which are then used to determine the 
corrections for the function 7 introduced in Ref. Eqs. (17), (30). (Also see below for 
definitions.) The first order terms in the S-matrix expansion for the self energy vanish, and 
so we consider the second order direct and exchange diagrams (see Fig. 4). These terms are 
dominant due to the heavy mass of the cantilever^. A three dimensional calculation for the 
Green's function is used to provide a more realistic situation than the one dimensional model 
introduced by Caroli et al. Due to the translational symmetry in the directions parallel to 
the barrier, we Fourier transform our Green's functions in these directions, and we obtain 
the expressions for /(V) by integrating over the distribution of possible fey-vectors. For the 
dispersion we assume a single parabolic band and leave the electronic rest mass unaltered 
for simplicity. We adopt the notation used by CCNS, that is: 

G — total propagator of the coupled system including phonon corrections; 

G° — propagator of the total system without phonon corrections; 

g — propagator of the right hand side infinite medium including phonon corrections; 

g0 — propagator of the right hand side infinite medium without phonon corrections; 

g — propagator of the left hand side infinite medium including phonon corrections; 

g° — propagator of the left hand side infinite medium without phonon corrections; and 

7(^0) = hm^'^o d x d x >g(x, x'\ k\\;u) (with 7 defined accordingly). 

The quantity 7 is related to the Green's functions of the semi-infinite media, which was 
introduced in Eq. (17) of Ref. |2|; we refer the reader there for a detailed treatment. 

The corresponding expressions for g° and g° were derived in the appendix of Ref. [| In 
our case the three dimensional treatment alters the u dependence slightly. The resulting 
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differential equation for the semi-infinite medium on the right hand side, after Fourier- 
transformation in the k» direction, now reads 

_ ^ + Vr{x) ~ {uj+ ^) 9 ° r{x > X] h]Uj) = 6{x ~ x) (3) 

where the superscript r denotes the retarded function and V r (x) has the shape shown in Fig. 
3. 

Explicitly, the unperturbed retarded Green's function g 0r is given by 

sinh K x 

—2m — — — — [(K + iq) exp{K(x' — a)} — (K — iq) exp{K(a — x')}] x < x' < a 
KD 
sinh K x 

—2m — J^Y) — ^ exp{K(x — a)} — (K — iq) exp{A"(a — x)}} x' < x < a 

sinh Kx' , . . . , 

—2m — — — — 2K exp{tq(x — a)\ x>a 
KD 

(4) 



g 0r (x,x'; k\\;u) 



k 2 / fc 2 

where K = \/2m(Vo — {u + // — ^)) and g = v/2m(a; + /i — ^) denote the wave vec- 
tors inside and outside the barrier respectively, and the denominator D is given by 
.D = (K + ig)exp(— ifa) + (K — ig)exp(A"a) = 2(K cosh(Ka) — sinh (if a)). For small 
barrier transparencies, which are considered here, exp(Ka) ^> 1 and the denominator is 
then approximately D = (K — iq) exp Ka. 



III. SELF ENERGY CORRECTIONS 

A. Dyson equation for the self energy 

Since our system is not translationally invariant we must resort to a Dyson equation 
in coordinate representation for the direction perpendicular to the barrier . We make the 
approximation, which will be justified in the following, that only the electronic subsystem 
on the right will receive any corrections from coupling to the mechanical degree of freedom. 
As the calculations for the unperturbed Green's functions on both sides are analogous, we 
will confine ourselves to calculating g, the propagator of the right subsystem. Ultimately, 
as shown in Ref. |2|, it is necessary to obtain corrections to the Green's functions g 0+ and 
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g° occurring in the Keldysh formalism in order to calculate a correction to the current. 
They are defined in Keldysh§ and also in Eq. (15) of Ref. ^| as 

g 0+ (u) = 2tW)/V) (5) 
g°-(u) = 2irip (u){f(u) - 1} . 

Here p°(w) and f°(oj) are the spectral density and the Fermi occupation function of the right 
subsystem, respectively. However, in the quasi equilibrium situation considered here, these 
can be obtained from the renormalized retarded Green's function g r . The corresponding 
Dyson equation reads 

g r {x,x';k { \-u) = g 0r {x,x';kfu) + J dx"dx'"g 0r {x,x"; k f u)t r (x", x"'-k\\\ u)g r (x"',x'; kf u) 

(6) 

Due to the localization of the electron-cantilever interaction the self energy has the form 

YT{x" ', x'"; k\\] oS) = 5(x" — a)5(x'" — a)S r (a, a; fcn; to) (7) 
and the integration simplifies to 

g r (x, x'\ fey; uj) = g 0r (x, x'\ k\\; oj) + g 0r (x, a; fcy; a;)E r (a, a; k\\;w)g r (a, x'\ k\\] uj) . (8) 
Expressed in terms of the unperturbed Green's function this is 

1 — g (a, a; kf u>)n r {a, a; kf oj) 

B. The self energy in Migdal's approximation 

In our system, the cantilever mass m c is much larger than the electron mass m e so 
that Migdal's approximatiori0 holds well. The self energy corrections will be small and the 
largest contribution to the self energy will come from the second order direct and exchange 
diagrams (see Fig. 4). Higher order diagrams will only contribute as , which in our case 



is certainly small. In view of the fact that we are calculating the lowest relevant orders of 
the corrections, we also neglect the second term in the denominator of In the Matsubara 
representation the direct and exchange diagrams are 



(2) . (ftk' 

£ ld (a, a; fc„; iu n ) = -(x V ) 2 D°(0)— J ^-jL £ G°(a, a; k[-iu' n ) (10) 



and 



E lex (a, a; k\\]iuj n ) = -(x V ) 2 ^-^2D°(iuj' n )G (a,a; k\\;iuj n - iu'J . (11) 

In ( |IUD T} d does not depend on k», i/W is the cross sectional area of the junction and 
/3 = (fcsT) -1 . The single mode phonon propagator is defined as 

The fact that there is no integration over the parallel momentum in the exchange term 
results from the fact that the electron-cantilever interaction is invariant in the transverse 
direction. The direct term from ( JTO| ) is not of any further interest, since it just renormalizes 
the equilibrium position of the cantilever, and we assume that this correction is already 
included in the definition of the distance 2a. Due to the symmetry in the problem the wave 
function \% of the entire system including the rigid square barrier factorizes as 

Xi(r) = <f>k x {%)(k l{ (p) (13) 

where p labels the coordinates in the parallel directions in real space and (p) is just a two 
dimensional plane wave. We can use this to perform the frequency summation in ( |TT|) using 
the representation for G° (cf. Eq. (28) of Ref. |2|), 

One further conceptual difficulty in calculating an exchange interaction diagram is that an 
electron can, in principle, travel from the position of the cantilever through the barrier to 

8 



the fixed electrode, and back again. This is why we have to take the propagators of the 
entire system as indicated in (|i0|) and (0). There is an analogous correction to the left 
side propagator g°, since an electron could travel to the oscillator position on the right 
and back again. These processes would allow the full semi-infinite propagators g and g to 
pick up a dependence on the chemical potential on the left hand side and right hand side, 
respectively. For the small barrier transparencies considered here, contributions arising from 
these processes can be shown to be negligible (see also sec. 3 of Ref. |]). To a very good 
approximation we can therefore replace G by g in (|10D and (0). We also replace in (|14|). 



the k -sum representation for g, the full wavefunctions ip by the wave functions 0, which are 
the analogous solutions for the semi-infinite problem. These wavefunctions are given by 

sinh(i^a;)e~^ a < x < a 

H D i5x a) 2\R cosh Ka + jq sinh Ka) ^ 



D 



a < x 



where now K = y2mVo — k 2 , q =\ k x | and D = 2 (K cosh Ka — iq sinh K a). The k -sum 
representation for g is then 



4klsmh. 2 (J2mV - k 2 x a) \ 

g (a, a; k»; iu n ) — 2 ; = • (16) 

Z 2™V -k 2 x + 2mV sinh 2 ( J2mV - k 2 x a) - H 



In order to simplify the notation, we omit the dependence of our functions on k\\ and a for 
the moment, and we set 



Akl smh 2 U2mV -k 2 x a) 

f(k x ) = ^ . . (17) 

2mV -k 2 x + 2mV smh 2 (J2mV - k 2 x a) 



(18) 



The Matsubara summation in (|TTJ) is standard and yields 
X 1 ™(iu n ) = -(x Vo) 2 2 Zf(k 



riB^c) + n F (£k) n B (u c ) + 1 - n F (e k 



iuj n + u c - ej: iuj n - u c - e k 



where = ^- — \x and np and n B are the Fermi and Bose distribution functions 



i 



If (JT8]) is considered at zero temperature the Fermi distribution functions turn into 6- 



functions and the Bose contributions vanish, so that 



9{-e % ) 9(e- h 



(19) 



As it turns out, the zero temperature approximation is not very accurate for realistic pa- 
rameters of a model system. The effects of residual temperature will be considered more 
thoroughly in the discussion and the appendix. 

The result of fliTf) can be analytically continued ( iu) n — > u + i5 ) to yield the retarded 



function. 

It should be remarked at this point that within the full Keldysh formalism for fully 
non-equilibrium situations, the Dyson equation (f|) in its correct form would read 

g + = (1 + + 9 a ^ a ) + <7 r £V (20) 

where the superscripts r and a refer to the advanced and retarded functions respectively (cf. 
Eq. (5a) in Ref. |3|). Also the retarded self-energy would not just be a convolution of G r and 
D r , but (cf. Eq. (6) in Ref. g) 

S r oc (D r * G-) + (D+ * G r ) . (21) 

However, in our quasi-equilibrium case both equations return to their usual equilibrium 
form. The main point in which the CCNS formalism enters our treatment is that we will 
renormalize the quantity 7 defined above to obtain a correction for the current in the next 
section. It is well known that the tunneling current across the unperturbed barrier is a non- 
linear function of the applied bias V for large enough biases, although the expression can be 
linearized for small biases. The modulation of the barrier width introduces a second, small 
energy scale, so that the resulting contribution to the differential conductivity significantly 
depends on V even for small values of the bias. 

Following an approach by Rickayzen0 and Scalapino^l the k -sum is assumed to have 
its strongest contributions coming from the immediate vicinity of the Fermi surface. Going 
over to an integral representation of the sum and changing the sign of in the first term 
of (U), we find 
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ilex 



(«) = -(x V ) 2 iV 1D (0)/(^/2m /U -^) ^ 



1 



+ 



1 



u + u> c + e + iS to — tu c — e + iS 



de (22) 



where iV 1D (0) 



2vr 1 



fi — k-IL) 1 is the one dimensional density of states at the Fermi 



surface including the sum over two spin directions. 

At this point it would be natural to attempt to include the effects of finite mechanical 
damping in the cantilever motion. However, it is difficult to include such effects from first 
principles and beyond the scope of this first approach prescribed herein. Instead, a simple 
way to model the effects of damping is to change uj c to the frequency for a damped classical 
harmonic oscillator, satisfying the equation of motion 



x + —x + uj n x = 

Q 



(23) 



which has the solutions u — ^ ± cu c y 1 — This replacement just introduces a finite 
imaginary part into the Green's function to achieve line-broadening and finite peak heights 
related to the Q of a classical oscillator. For simplicity let b = uj c ^Jl 
The self energy can then be written 



W and c " ~ 



2Q • 



ilex , 



u) = (x V ) 2 NfZ\u;) 



(24) 



where 



de 



(oj + b + e) — ic 
(u + b + e) 2 + c 2 



IC 



(u-b 



e) 2 + c 2 



b„ 

2 



(lo - bf + c 2 



(u + b) 2 + c 2 



7i — arctan(- 



arctan(— - — ) 



(25) 



The real part of S 6 (cu) in (25) has the usual logarithmic form for this kind of diagram and 
the imaginary part is a smoothed-out version of the 9 -functions occurring for zero damping. 



IV. CORRECTIONS FOR THE CURRENT 
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A. Derivation of the corrections 



The self energy obtained in the last subsection can now be inserted into the simplified 
version of ([|) to yield a correction for the right hand side propagator g, which is needed to 
calculate the corrections to the current and its first and second derivatives. We will use the 
same representation of the current as in Eq. (41) of Ref. ||] with the only changes due to the 
three dimensional nature of our model, so that the expression for the total specular current 
I x across the barrier is proportional to the cross sectional area of the system. Dropping the 
subscript, it reads 

wo> - ^ /: i / St-'t^B ms - - ^ - «■» < 26) 

where uj = + /iR = uj^ + fi^. For T = this can be recast into 



(27) 



In the limit \ll — > — > /i (^7|) goes to 



(/(V)) - 8 ^ /„ 2^7 (2tt) 2 I 7 (0) + 7(0) | 2 (28) 
(see the discussion in Ref. ||] for the derivation). In the above expressions the imaginary 
parts of 7 still contain the k\\ dependence, which is not indicated for simplicity. According 
to its definition above, 7 on the right hand side will receive a correction through the Dyson 
equation for g, yielding 7 = 7 + 7' , whereas for the left hand side there is no correction, so 
that 7 = 7 . Following CCNS it is assumed that for small bias 7 ~ 7 holds well. In order 
to get the first order correction to the integrand in fl28|) it is expanded to first order in 7', 
yielding 

Jm7°(/m7 + Imrf) (Im-/ ) 2 | Imj°Imj' 4(Im'y ) 2 (Rej Rej' + Imr^Iwrf) 
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I 7 + ( 7 + y) |2 J 2 7 |2 I 27O |2 ' J 27 

(29) 

The expression for 7 is already given in Ref. § and can be written 
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iq cosh Ka — K sinh Ka — (q 2 + K 2 ) cosh Ka smb. Ka + iqK 
if cosh if a — iq sinh if a X 2 cosh 2 Ka + g 2 sinh 2 Ka 

For arguments ifa ~ 4 or larger, which seems reasonable for real tunneling processes, ( [301) 
goes asymptotically to the WKB-like result 



7° ~ 



-K + 



AiqK 2 
2mV ' 



(31) 



This shows how the exponential dependence of the tunneling current on the barrier width 
enters this calculation. 

The correction 7' is then obtained from the definition of 7 by differentiating the correction 
term in the Dyson equation 
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2m x,x'->o 



lim d x d x ig°(x, a)S(a, a)g°(a, x') 

>y ! — *n 



8mK 2 
D 2 



E(a, a) 



(32) 



where again the dependence of g on fen is not indicated explicitly. For reasons that will 
become apparent as we proceed, it is convenient to calculate the dimensionless ratios 
and . Written in the following compact way we obtain 

Irwy Rey 



-x z Q V (2m) 2 lV o -(u + fi-^)\^0i 



2m' 



Im(£) . i?e( s 



D 2 - 



7m7° _Re7° 



(33) 



Using the expressions introduced for 7 and -D , fl33|) then yields 



7m7' J i?e7' 
irrry Rej 



C 



K 2 - a 1 

ReE b + ImE b 



2qK 



2 {2qKImL h - (K 2 - q 2 )ReT} 



-2Ka 



2mV 



(34) 



k 2 



where C = ^mx^ ^/i — ^ J ^Vo — (a; + /x — • Equation fl34|) clearly shows that the 

second part of the square bracket is negligible for the small barrier transparencies considered 
here. 

We can now conveniently write 
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(35) 



where we define rj = j^o ■ Inserting (|35|) into ( f29|) now shows that the third term in fl29|) 
arising from the expansion of the denominator will be exponentially smaller than the first 
two, due to the exponential dependence of Irwy . Thus (|29|) can be written 

Jm7°(Jm7° + Irnrf) (Jm7°) 2 . 
| 7 o + ( 7 o + y) [2 - | 2 7 o p ( +r l) ■ 

So far all our results depend on k\\, and as CCNS mention, the unperturbed expression 
in the integrand of ( p8|) is a strongly peaked function around ku = 0, which mainly comes 
from its exponential dependence on Ka . More precisely, differentiating with respect to V, 
(m yields 



di\ r<Pk\ U ,„ fcf, 

dvl=i J^f h{eV + ^-2^^ V \- Aa \ 



( k 2 
2m iv Q - {eV + fi - JL) 



(36) 



where h is a function that varies only slowly over the range of the integration compared 
to the exponential. Due to the strong exponential behavior of the final term, the largest 
contribution to this integration will come from the vicinity of the lower integration limit 
after transforming to polar coordinates. The square root in the exponential in (|36"D can be 
expanded in z = k 2 around ku = 0. The integrand can then be cast into: 



' dV 

dV. 



00 dz . 



h(eV + a ) exp <^ -AaJ2m \V - (eV + u)] > exp 

« Itt 2m ] P l V m \ I . /2m [V - (eV + //)] 



-2az 



pm [V Q -{eV + n)] 
8iia 



h(eV + jj) exp <^ -Aad2m [V - (eV + fi)] 



where we have neglected the contribution coming from the upper limit. 



B. Discussion of I(V), dl/dV and d 2 I/dV 2 

The current- voltage characteristic /(V) can now be obtained by integrating (|36|) . For 
two reasons this procedure is not very rewarding. Firstly, the integral is difficult to perform 
as the exchange term is no longer a slowly varying function over the range of integration, 
and thus our previous approximation of replacing J^ +Sx f(x')dx' by f(x)5x is not accurate. 
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Secondly, the interesting corrections to I(V), which come from the exchange diagram are 
already very small in dl/dV (as will be shown later), so that the integration leading to 
I(V) will make them essentially undetectable. We will thus concentrate on the expressions 
for the first and second derivatives of /(V). The first derivative dl/dV, i.e. the differential 
conductivity, is simply given by Q3"7|) and accordingly d 2 I/dV 2 can be obtained by a further 
differentiation. We thus have 

eV 2 ) 



dl 

dV 



2ji 2 a 



(7m7 ) 2 



27 ! 



d 2 I 



;V 2 ) d 



(38) 



dV 2 2n 2 a dV 

As an approximation to (im 7 ) 2 / \ 27 \ 2 we have, consistent with our approximation in 
(|36|), taken over the WKB-expression (for k\\ = 0) derived in Eq. (47) of Ref. for further 



calculations, 

(im 7 ) 2 



(eV + n) [Vo-jeV + fj)] 
V 2 



exp |-4 av /2m [V - {eV + //)]} . 



(39) 



| 27 | 2 

We have used the expressions derived above to calculate the magnitude of these effects 
on the I-V characteristics of a model system, consisting of a single-crystal Si beam, rigidly 
clamped at both ends. The beam is chosen to have a length of 500A and a cross-section of 
100A x 100A, giving it a fundamental resonance frequency of uj c /2tt = 30GHz and a mass of 
1.2 x 10~ 20 kg. The effective spring constant of the beam will be 420 N/m. The parameters 
for the calculation are listed in Table 1. 



Barrier 


Vo = 5eV 


fi = 2eV 


a = 5A 




Cantilever 


u c = 1.2 x 10~ 4 eV 


m c = 1.2 x 10~ 20 kg 


k c = 420iV/m Q = 10 3 



It should be noted that these results are calculated at T = 0, and that a rather modest 
quality factor of 10 3 for the oscillator has been chosen. Mechanical oscillators at much lower 
u) c , but with Q's of order 10 6 or higher have been reported^. Our choice of small Q was 
made to very roughly compensate for the finite temperature broadening effects in any actual 
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experiment. As described in the Appendix, a Q of 10 3 corresponds roughly to a measurement 
on this system at lmK, but with a mechanical Q possibly much higher than 10 3 . 

Figures [5] and |6| show the first and second derivative of the current density j (in A/m 2 ). 
The unperturbed terms for dl/dV and d 2 I/dV 2 are not symmetric about the origin, since in 
our model we have assumed that the Fermi energy is fixed on one side of the barrier, while 
it is varied about this value by the bias on the other side. It can be seen that the relative 
magnitude of the exchange contribution in dl/dV to the unperturbed value is about 1CT' 
if temperature effects are neglected entirely. The form of the exchange part is similar to 
that observed for localized phonons. Comparing these graphs, it is apparent that it would 
be difficult to see the contribution of the exchange resonance after a final integration to 
get I(V). The most visible feature is the correction to d 2 I/dV 2 : the exchange correction 
here reaches a size comparable to the unperturbed value, and the resonant form clearly 
distinguishes it from the smooth background. 

Since the graphs only give a qualitative impression of the peak strengths, we now investi- 
gate them more thoroughly. The peaks in dl/dV clearly result only from contributions of the 



real part of the exchange term, cf. Eq. (^Sj). The term has poles at uo = ±b = ±w c wl — 
which, in the bare term, leads to a peak height of i?eS b (±6) = =p| [In 16 + InQ 2 ], indepen- 
dent of u c . The origin of the peaks in d 2 I/dV 2 is slightly more delicate. Here, the largest 



contributions come from the terms including -jj- 
term has poles at 



ReZHu, 



and also -r- 



. The first 



u = ±\lb 2 + c 2 ±2cVb 2 + c 2 = ±ujJ 1 ± — ~ ±u c (l ± — ) , 

where the choices of the signs are made independently. The peak heights scale as ±Q/u c 
for the outer (inner) poles respectively. The second term has its poles exactly at u — ±6 , 
with heights that scale as ±2Q/uj c , and width (FWHM) of u c /Q, and thus the peaks in the 
second derivative do depend on lu c . To see how Q quantitatively affects the peak heights, 
we have compiled the following table of maxima of the peaks in absolute height occurring in 
the exchange corrections to dl/dV and d 2 I/dV 2 , where apart from Q all other parameters 
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for the system remain unchanged. 



Q 


10 3 


10 4 


10 5 


10 6 


dj/dv[n- l m- 2 } 


150 


188 


227 


265 


d 2 j/dV 2 [G(QV)- 1 m- 2 ] 


0.23 


2.1 


21 


210 



Eventually, as the intrinsic Q of the oscillator is increased, the broadening at zero tem- 
perature will be determined by electronic damping of the oscillator, given by higher order 
terms in the self energy. However for the values of Q in the table, we believe these effects 
are unimportant. 

V. CONCLUSION 

In conclusion we have investigated the first and most relevant self energy corrections to 
the first and second derivatives of the I(V) characteristic in M-I-M tunneling due to the 
interaction with a single localized mechanical mode of a movable tunnel junction. This 
represents the first step in investigating the influence of, and back-action on, oscillatory 
mechanical degrees of freedom in a mesoscopic tunneling system. We find that the presence 
of this mechanical mode gives a pronounced signature in dl/dV, and a much stronger one in 
d 2 I/dV 2 , in the regions where the bias across the junction equals the energy of the eigenmode 
of the oscillator. This strong non-linear enhancement of the differential conductivity and its 
derivative across the junction in this region can be viewed as arising from the opening of a 
new phonon-assisted channel for inelastic tunnelling. 

As the derivation of these results has been made using the general many body formalism 
of CCNS, the approach may be readily extended. Possible applications might include con- 
sidering stronger bias to generate a fully non-equilibrium situation, incorporating the band 
structure in the metals, and elucidating the role of other many body interactions. 
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VII. APPENDIX 



In this appendix we derive the range of validity of the zero temperature calculation 
and establish qualitatively the changes to our results as finite temperature effects start to 
compete with and finally dominate the intrinsic losses of the mechanical oscillator. The first 
zero temperature approximation which entered our calculation from before was in (|i~8|) for 
the Matsubara self energy. 



-ilex / 



lUJr 



n B {oJ c ) + n F (e^) n B (u c ) + 1 - n F (e- k 



iu n + U) c - £g lU) n — u c — eg 



(40) 



If one still assumes that the largest contributions to the ^-suni come from the vicinity of 
the Fermi surface one can approximate the sum similar to (E3) as 



Ej"» = -(x V ) 2 N 1D (0)f^2m^~k 2 ) 



x ( n B (co c ) I de 

-oo 



1 



+ 



lu + lu c — e + i6 lu — uj c — e + iS 



- I de 

— oo 



n F (e 



+ 



1 — npfe) 



lu + uj c — e + i8 uj — uj c — e + iS 



The two integrals can be considered separately and the first one gives 



de 



1 



1 



P 



de- 



2(oo-e) 



to — e) 2 — oo 2 



2m 



(41) 



(42) 



_uj + ijj c — e + iS uj — <jj c — e + iS . 
where p denotes the principle value. The principle value integral is seen to be zero, since 
the integral is antisymmetric. The imaginary part turns out to be constant with respect to 
u. The second integral can be recast into 



de 



2lu 



in I de 

If that is compared to formula (0), one sees that the difference for finite T is that the 
integrals now get a smooth cut-off at e = as opposed to the ^-function cut-off from before. 
On introducing a finite quality factor Q, fl43|) goes over to 

de 



8(e + lu + uj c ) + 5(e — lu + lu c 



(43) 



e-P £ + 1 



(lu + b + e) — ic (u> — b — e) — ic 



(44) 



(lu + b + e) 2 + c 2 (u-b-e) 2 + c 2 _ 

in analogy to (|25|) . We can thus choose either ( ^3l) or (0) as a basis of reasoning, depending 
on the regime considered. As it turns out, for the conceivable range of oscillator quality 
factors Q, the linewidths induced by Q will be very small compared to the temperature 
broadening, even for temperatures as small as ~ lmK. Thus, we assume for now that 



Q — ► oo and use formula ([43"D. The difference in the behavior for the real part of (f43D above 
from its T = analog can be seen through partial integration: 

1 2lu p r°° de 



4 ./-oo cosh 2 f £ 



111 



LU — 10 r 



lu + lu c + e 



-oo e-^ + 1 cj 2 - (lu c + e) 2 
where the integrated term has vanished at the boundaries. In the limit (3 —>■ +oo 

(3 1 



(45) 



4 cosh 2 f £ 



5(e) 



(46) 



This way one can see how the T = result for the real part of (51) is recovered. The 
corresponding integration of the imaginary part gives: 

1 1 



ITi 



+ 



e /3(u>+a; c ) _|_ J g— f3{u—u c ) _)_ ]_ 



(47) 



which is seen to have a smoothed out step behavior, being essentially constant for | lu |^> cj c 
and being more or less zero for | lu |<^C lu c as is seen in the corresponding expression in 
(Pl|), with the only difference that the width of the transition is now the width of the Fermi 
distribution, which is approximately 
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3.5k B T ~ 3 x 10~ 4 T eV/K . (48) 

It is possible to determine phenomenologically what effect this will have for the most "ob- 
servable" of the peaks, i.e. the ones in d 2 I/dV 2 . The peaks in that function are, as we 
found, determined by dY} r ex (u) j du . It can be seen directly from (|47D that the derivative 
of the imaginary part, which gives the strongest contributions to these features, exhibits a 
pair of Lorentzian peaks at the resonance frequency of ±a> c . At nonzero temperature these 
have a width determined by the Fermi distribution, as shown in fl48|) . One can also see the 
change in the behavior of the derivative of the real part of E from fl45|) without having to 
perform the integral. Differentiating with respect to u gives 

4 7-oo cosh 2 § eoo 2 - (u c + ef ' 1 ' 

This can be viewed as being a sum of two convolutions of opposite sign in u between the 
derivative of the Fermi function 1 — rip and the two partial fractions of the second factor. 
The result of the convolution of these two peaked functions is to a good approximation again 
a peaked function, where the width of the resulting feature is the sum of the widths of the 
convoluted components. Since for the moment we assume that the width of our initial peaks 
(due to the quality factor Q) is practically negligible compared to the thermal width of the 
Fermi function, the effective width of the peaks will also be given by the thermal width. If, 
on the other hand, in an intermediate regime (low Q, low T) the two effects start to compete 
one would have to use flUf) for a correct treatment. However, as mentioned above, the effect 



will just be that the widths will add to give the result. 

The second T = approximation was introduced in ( pTf) for the current. In analogy 
to the reasoning above, one can argue that if the Fermi factors are kept the smoothed out 
cut-off will lead to similar effects, which will enhance the total broadening by an additional 
factor of about two. In summary, we can conclude that the contributions of the quality 
factor Q and the temperature T to the total linewidth will be about 



8uj 



land 



8uj 



U)r. 
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= 6X1^, 



and the total width will be approximately the sum of these two contributions. 
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FIGURES 

FIG. 1. Schematic view of tunneling electrode on cantilever with mass m c and spring constant 
k c , placed a distance 2a from an infinitely massive counterelectrode. The device is biased with a 
voltage V and the resulting current I is measured. 

FIG. 2. Full potential barrier used in the calculation. The tunnel barrier has height Vq and 
equilibrium width 2a. The right hand side of the barrier corresponds to the position of the cantilever 
tip along the x-axis. 

FIG. 3. Semi-infinite potential barrier used to split the Hamiltonian into two pieces. The 
potential goes to infinity at x = 0, is equal to Vq for < x < a, and is zero elsewhere. 

FIG. 4. Lowest order diagrams contributing to the self energy: (a) direct; and (b) exchange. 
Electron paths are shown by straight lines and cantilever "phonon" paths are shown by sinuous 
lines. 

FIG. 5. Correction to the first derivative dJ/dV of the current density: the solid and dashed 
line are for cantilever quality factors Q = 10 3 and 10 respectively. 

FIG. 6. Second derivative d 2 J/dV 2 of the current density containing the corrections considered 
in the text; the solid and dashed lines are for Q = 10 3 and 10 respectively. Inset: second derivative 
for Q = 10 3 , with vertical scale sufficient to show full extent of peaks. (Axes have same units as 
main plot). 
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